Concentration of Charge Carriers and Anomalous Gap Parameter in the Normal 

State of High-T c Superconductors 
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Fermi-Dirac statistics has been utilized by introducing the 
average ionization energy (Ei) as an additional anomalous 
energy gap in order to derive the two-dimensional concentra- 
tion of charge carriers and the phenomenological resistivity 
model for the superconducting polycrystalline materials. The 
best fitted values of Ei and the charge carriers' concentration 
ranges in the vicinity of 4 to 9 meV and 10 16 m -2 respectively 
for the superconducting single crystal samples and polycrys- 
talline compounds synthesized with various compositions via 
solid-state reactions. The phenomenological resistivity model 
is further redefined here based on the gapless nature of charge- 
carriers' dynamics within the CU-O2 planes that corresponds 
to anomalous Fermi liquid behavior, which is in accordance 
with the nested Fermi liquid theory. 

PACS Number(s): 74.72.-h; 74.72.Bk; 71.10.Ay; 72.60.+g 



I. INTRODUCTION 

Since the reported discovery of high-T c superconduc- 
tors (HTSC) based on the copper oxide materials, the 
scientific communities have been swept with revolution- 
ary excitement. This leads to intensified focus on this 
particular material both in term of theoretical and ex- 
perimental research throughout the world. The discov- 
ered Ba-La-Cu-0 system by Bednorz and Muller Q with 
a T c of 33 K have had led to subsequent findings of su- 
perconductivity in La-Sr-Cu-0 at 40 K (||| an d later 
in Y-Ba-Cu-0 system at 90 K [@. Actually the La- 
Ca,Sr,Ba-Cu system were first synthesized and studied by 
Raveau et al. H • The conduction mechanisms in the nor- 
mal state of HTSC such as the out-of-plane and in-plane 
resistivities are intriguing due to the highly anisotropic 
behavior. Anisotropic normal state electrical properties 
such as metallic-like conduction in the CU-O2 planes and 
semiconducting-like conduction between those planes || 
are believed to be one of the central concerns for ap- 
plications 0. This in turn may lead one to extract 
some clues to understand the mechanism of high-T c su- 
perconductivity Interestingly, Levin et al. |J us- 
ing the T-dependence on anisotropy p c /p a b showed the 
conductivity relationship between cr c and a a b as: cr a b — 



(a + bT~ 2 / 3 )a c . This T~ 2 / 3 dependence point towards 
a 2-dimensional (2D) localization and hopping. Quite 
recently, Basov et al. |To| have analyzed the interlayer 
infrared conductivity of HTSC. They found that discrep- 
ancy occurred between superfluid density, p s and differ- 
ence between normal state (N n ) and superconducting 
state (N s ) spectral weights i.e. p s > N n — N s . This 
discrepancy was only revealed for c-axis electrodynamics 
but could not find a similar conflict in a6-plane (p s = N n 
— N s ), which can be interpreted as a gapless behavior 
of charge carriers in a6-planes. PrBa2Cu307 compound 
gives a unique non-superconducting, semiconducting-like 
normal state electrical property and the nature of this 
unique property was discussed intensively in terms of 
hole doping in CU-O2 planes [p~T|— p~3|| . As a matter of 
fact, semiconducting-like behavior of p(T) corresponds 
to localization of charge carriers connected with oxy- 
gen disorder in G\x-O x chains and the disorder in Pr 3+ , 
Pr 4+ distribution fl4f| . Furthermore, positron annihila- 
tion studies Jl5|,[l6| showed that PrBa2Cu3C>7 compound 
has the Fermi surface due to Cu-O^ chains' state. The 
concept of charge concentration in high-T c cuprates is 
very ambiguous [pL7|| . However, it is well established that 
T c of cuprates depends on holes' concentration within 
CU-O2 planes where carrier concentration is determined 
by counting the valence number frij , |l8|| . Very recently, 
concentration of interlayer electrons and its resistivity 
model for HTSC was derived via Fermi-Dirac statistics 
with an additional gap constraint, which is the average 
ionization energy (Ej) jl^. In this letter however, a phe- 
nomenological polycrystalline resistivity model is derived 
based on E] as well and is redefined here to obtain an 
a6-plane resistivity model by incorporating gapless na- 
ture of charge carriers in CU-O2 planes. Besides that, es- 
timation for concentration of charge carriers and Ej via 
best fitting in the normal state of HTSC are also illus- 
trated in this work. Apart from that, estimation of resis- 
tivity transition such as from metallic to semiconductor 
or vice versa upon substitution in crystalline 123 phase 
Pr a; Yi_ 2: Ba2Cu307 system is also highlighted. Ions' av- 



erage ionization energies were calculated from ref. 20 1. 
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II. CHARGE CARRIERS' CONCENTRATION 
AND RESISTIVITY MODELS 

Tunneling electron's distribution has been derived via 
Fermi-Dirac (FD) statistics using average ionization en- 
ergy (Ex) as an anomalous constraint, and is given below 



f e (E) = exp[-n-\(E + E I )}. 



(1) 



Similarly, probability function for holes can be derived 
and is given in Eq. |^ below, where it is obtained by in- 
serting the restrictive condition E hole = E mitia i state - 
Ei. Einiuaistate is the energy of a particle in a given sys- 
tem at a certain initial state and ranges from to -oo. 
This can be justified as follows; for a hole to occupy a 
lower state M from initial state N is more probable than 
to occupy state L if the condition Ej(M) < E](L) at 
certain temperature, T is satisfied. 



f h (E) = exp[n + \(E - E T )]. 



(2) 



Note that Eq. [j] and Eq. || can be reduced to the stan- 
dard Fermi-Dirac distributions if Ej approaches 0. The 
importance of average E^s inclusion is that, it can be in- 
terpreted as a c-axis gap, which represents the gap that 
tied to a neutral atom and could be varied with ionic 
substitutions. The respective normal state concentration 
of electrons [Ell and holes are given by 



P 



f e (E)N e {E)dE, 



f h (E)N h (E)dE. 



(3) 



(4) 



fi = —Ep/ksT, A = 1/fcgT are carefully chosen to 
confine the total number of particles in a given system 
equals to n. N e (E) and Nh(E) are 2D density of states 
(DOS) that can be derived from the Schrodinger equa- 
tion. N e (E) and Nh(E) are derived as ml/nh 2 and 
m* h /ith 2 respectively and note that both 2D DOS are 
independent of E. Solutions of integral Eq. || and Eq. ^ 
give the concentration of tunneling electrons between Cu- 
O2 layers, n and concentration of holes, p in the normal 
state of HTSC as 



n = = — exp 

irh 2 



m* h k B T 
irh 2 



-exp 



Ep — Ei 
k B T 



-E F - Ei 



(5) 



(6) 



Eq. |^ and Eq. ^| are expressed in term of ionization 
energy, where fee is Boltzmann constant, m* and m* h are 
electrons' and holes' effective masses respectively, h = 



h/2%, h is Plank constant and Ep is Fermi level. More- 
over, from Eq. ^ and Eq. ^, concentration of charge car- 
riers in the normal state of polycrystalline HTSC can be 
written in the form of geometric mean which is given by 



np 



k B T 
nh 2 



yjm* e rn* h exp 



-Ej 
kuT 



(7) 



This equation is free from Fermi level, hence concen- 
tration of charge carriers in superconducting polycrys- 
tals can be calculated by first estimating Ei, which is 
described in the following paragraph. Scattering among 
fermions leads to the conclusion of 1/r oc T 2 M and has 



been confirmed empirically that p oc T as well [|21|-[23|. 
Furthermore, by utilizing metallic resistivity equation of 
p = m/ne 2 T where e and t are electron's charge and mean 
scattering free time respectively |24|] then one can derive 
the normal state phenomenological model for HTSC as 
given in Eq. ^ below. This were carried out by substi- 
tuting m — to* = m* h , n — y/np, and 1/r = AT 2 into 
metallic resistivity equation. 



Ppoly 



A- 



-Texp 



Ej 



k R T 



(8) 



A is independent of temperature. Eq. |^ above can be 
used in the normal state resistivity of polycrystalline su- 
perconducting samples to evaluate the variations in Ej. 
Substituting this Ej into Eq. |?] at certain temperature, T 
leads to the estimation of charge carriers' concentration. 
Subsequently, an etfr-plane resistivity model can be de- 
rived quite naturally by taking anomalous gap = due 
to gapless behavior of charge carriers in CU-O2 planes 
ifiof , thus Eq. H can be simply rewritten as 



Pab (T) = A—T. 



(9) 



This latter phenomenon could be due to change in the 
scattering rate that results from a different Fermi sur- 
face phase space argument in accordance with the nested 
Fermi liquid theory |2q|. 



III. DISCUSSION 

Polycrystalline resistivities, p po iy{T) of crystalline 
YBa2Cu307 |2(| samples (B and C) were obtained us- 
ing 



Ppoly (T) = y/p c (T)p ab (T). 



(10) 



It is assumed that Eq. |To|is valid to convert p c (T) and 
Pab{T) of single crystals to p po i y (T) in order to analyze 
the effect or changes of A parameter between single crys- 
tals and polycrystals. Consequently, one can study the 
scattering effect solely on A parameter. Fitting of exper- 
imental Ppoly (T) with Eq. || gives respective values for 
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anomalous gap and concentration of charge carriers in 
the order of 4-9 meV and 10 16 m~ 2 as listed in Table Q. 
These fittings are depicted in Fig. [j] and Fig. || for sam- 
ples B and C respectively. Furthermore, resistivity curves 
of polycrystallineGal212 §7], Yi^Pr^SrBaCugC^ @, 
YBaaCuaCV, " 



and T11212 (3TJ samples are uti- 
lized by calculating the slope, dp(T)/dT in the normal 
state of respective samples. This slope would enable 
one to determine Ej via Eq. || and followed by deter- 
mination of concentration of charge carriers in the nor- 
mal state at 300 K and 150 K, which are also listed 
in Table 1. The above concentration of charge carri- 
ers was estimated by substituting the best fitted Ej into 
Eq. 0. All samples have charge carriers' concentration in 



the order of 10 m , however, parameters A and Ej 
vary accordingly with samples. Note that A is equals 
to AttH 2 /e 2 kB for simplicity's sake. It is quite interest- 
ing to note that all Y123 samples including the single 
crystals (sample A and B) have values for Ej between 
7.69 and 8.25 meV. Contrary to that, A parameter is 
approximately 100 times lower for single crystals (10 -8 ) 
than polycrystals (10~ 6 ) where this could be due to the 
nature of scattering in crystalline samples. Note also 



that the magnitude of p po i y (T) estimated from Eq. |lO|for 
single crystals are lower as mentioned previously. Simi- 
larly, A and Ej parameters for 1212 phase compounds, 
Tl(Sri. 6 Sm .4)CaCu 2 O7_ x , Tl(Sri. 6 Dy .4)CaCu 2 O7_ I 
H and GaSr 2 (Yo. 6 Cao.4)Cu 2 7 |7| are in the or- 
der of 10~ 5 and 5-6 meV respectively. Substitution 
effect in Yi-^Pr^SrBaCusOy system is found to be 
not in favor of superconductivity where, Pr substi- 
tution increases anomalous gap, Ej and A parame- 
ter. Parallel to this, comparison between polycrys- 
talline samples, YBa 2 Cu307_ x (Sample A) |^9| and 
Y(Ba .85Pbo.i5) 2 Cu307_2; (Sample B) Q expose a re- 
duction in Ej that suggests Pb substitution into Ba 
sites is somewhat favorable. Further comparisons be- 
tween samples prepared by different listed researchers 
are not discussed here since the exact synthesis method 
varies from one researcher to another. In another con- 
text, p c (T) transition that occurs in superconducting 
Yi_ x Pr x Ba 2 Cu307 system is in fact due to the differ- 
ences in average Ei of p r 3 +' 4 + and Y 3 +. The valence 
number, z+ (= 3+) and q+ (= 3+) correspond to ions, 
Pr z + and Y« + respectively. Since Ej of Y 3 + (1260 
kJmol" 1 ) > Ej of Pr 3 + (1210 kJmol" 1 ) then substitut- 
ing Y with Pr will reduce Ei — Ep which is the c-axis 
gap parameter [|l9| for single crystals thus, reduces p c (T). 
Similar substitutions would also decrease the Ei gap pa- 
rameter and Ppoiy (T) which is for polycrystals that can be 
verified from Eq. ||[ In contrast, the experimental data 
obtained by Jiang et al. (single crystals) |3l[] and Das 
et al. (polycrystals) [ fl3| expose a contradicting relation 
where p c (T) and p po i y (T) increases with Pr substitutions 
into Y sites. Thus valency of Pr is expected to be greater 
than 3+ due to Pr 4+ ions' contribution p2T|. Parallel to 



this, Meng et al. |3^] have employed PVL theory to de- 
termine valency of Pr in superconducting state. They 
have found that valency of Pr is vital in producing a su- 
perconducting Yi_ 2; Pr 2 ;Ba 2 Cu307 compound. In their 
calculations, valency of Pr is < 3.15 to get a supercon- 
ductor compound for all Pr doping but a valency of > 
3.15 gives a superconducting state for only up to a lim- 
ited Pr doping. Therefore, using the concept presented in 
section II, the minimum value of Pr's valence state, which 
is obviously above 3+ that contribute to non-metallic be- 
havior, can be calculated from linear algebraic equation 
given below 



- E u + - Y E u = - ^ E u 



(11) 



=2+1 



The first term, j £i=i+i E u in Eq. 



11 above has i = 
It is solely 



z + 1, z + 2,..., z + j and j — 1, 2, 3, 
due to Pr 4+ ions contribution or caused by reaction of 
the form Pr 3+ — electron — * Pr 4+ , hence j is equals to 1 
in this case and S represents the additional contribution 
from Pr 4+ . The second (i — 1, 2, 3, z) and last (i 
= 1, 2, 3, q) terms respectively are due to reaction 
of the form Pr — 3(electrons) — > Pr 3+ and Y — 3(elec- 
trons) — * Y 3+ . Recall that q — z — 3+ and i = 1, 2, 3,... 
represent the first, second, third, ... ionization energies 
while j = 1, 2, 3, ... represent the fourth, fifth, sixth, ... 
ionization energies. Therefore, z + S gives the minimum 
valence number for Pr which is calculated to be 3.0133 
from Eq. [ll]. This critical value corresponds specifically 
for p(T) variations in the normal state behavior where 
p(T) decreases with doping if and only if z + S < 3.0133 
and vice versa if z + 5 > 3.0133. Those changes in p{T) 
will eventually affect dp(T)/dT slope transition in term 
of metallic-like or semiconducting-like normal state prop- 
erties as mentioned for pure 1212 phase polycrystals [^9| . 
In order to justify this critical value, p c (T) measurements 
of Yi_ x Pra;Ba 2 Cu307 single crystals |3l[] that indicate a 
transition of dp c (T) / dT slope from positive to negative 
with x doping were used. Thus, from the above-presented 
approximation, z + S or valence state of Pr for x = 0.13 
to 0.42 samples is > 3.0133 since there is an increment of 
Pc(T) with Pr doping. In contrast, z + S value should be 
above 3.15 as mentioned previously in accordance with 
Meng et al. p3|| . Although 3.15 is larger than what have 
been anticipated (3.0133), but one should note that 3.15 
corresponds only to superconducting behavior without 
considering the normal state behavior. Therefore, by 
considering only the normal state behavior, decrement of 
p c (T) with increasing Pr doping only occurs if z + 5 < 
3.0133. In short, if 3.0133 < z + 5 < 3.15, then the mate- 
rials are superconductors for all x doping but dp c (T) / dT 
changes from positive to negative with x. However, if z 
+ 5 < 3.0133 then the materials are also superconduc- 
tors for all doping, x but dp c (T) / dT remains positive and 
p c (T) will decrease with x. Finally, ifz + £>3.15, then 
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the materials are superconductors for only upto a cer- 
tain x doping in which, dp c {T) / dT changes from positive 
to negative and p c (T) would increase with x. Note that 
the prediction above is only valid for very pure materials 
without any significant impurity phases. p c (T) data at 
high temperatures J3l[] surprisingly decreases with fur- 
ther Pr doping from x = 0.42 to 0.55 though it remains 
semiconducting. This phenomenon could be possibly due 
to variation in the magnitude of z + 8 where the value 
for z + S might change from; z + 5 (at x — 0.42) > z 
+ 5 (at x = 0.53) > z + S (at x = 0.55). I.e., Pr z+S 
varies with doping or concentration, x. Apart from that, 
as evidenced in experimental data, the phenomenologi- 
cal afr-plane resistivity model (Eq. ||) indicates a linear 
relationship between temperature and in-plane resistivity 
due to gapless nature of charge carriers, thus, pointing 
towards an anomalous Fermi liquid behavior within Cu- 
O2 planes. 

IV. CONCLUSIONS 

Two-dimensional Fermion characteristics have been 
employed to derive charge carriers' concentration. 
Anomalous gap parameter (Ej) and charge carriers con- 
centration that has been estimated via best fitting for 
all superconducting samples of polycrystalline and crys- 
talline are in the order of 4-9 meV and 10 16 m~ 2 respec- 
tively. As a matter of fact, variations that occurred in 
y / np, A and Ej are due to variations in sample's pu- 
rity, methods of preparation and its compositions. Phe- 
nomenological a6-plane resistivity model (Eq. ^) indi- 
cates a linear relationship between temperature and in- 
plane resistivity due to gapless nature of charge carri- 
ers, signifying an anomalous Fermi liquid behavior in 
accordance with nested Fermi liquid theory. In addi- 
tion, disorder of valence state in certain elements such 
as Pr 3+ , Pr 4+ in term of average Ej also contributes 
to the normal state behavior such as from metallic- 
like to semiconducting-like behaviors. In other words, 
dp c (T) I eff's varies from positive to negative or vice versa 
with selected doping or substitutions at certain temper- 
ature. Simply put, polycrystalline resistivity's slope is 
dependent on average ionization energy, where valence 
state of Pr which is above 3+ could be estimated from 
dp c (T) I dT transition upon substitution of Pr into Y sites 
as Eq. [ll] states. 
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TABLE I. Calculated values of the temperature indepen- 
dent parameter (^4), anomalous gap parameter (Ei), resistiv- 
ity slope (dp(T) /dT) and the concentration of charge carriers 
{y/rvp) at 300 K and 150 K in the normal state of high-T c su- 
perconducting compounds. Note that the magnitudes of Ei 
are in the range of 4-9 meV and the charge carriers' concen- 
tration is determined to be in the order of 10 16 via appro- 
priate fittings. dp(T)/dT is calculated from the experimental 
p P oiy(T) plots. Note that y/np is obtained from Eq. [?]. k B is 
the Boltzmann constant. 



FIG. 1. The resistivity curves of crystalline YBa2Cu307 
(Sample B: with volume of unit cell = 630 x 640 x 75 (im) 
of 123 phase were calculated using Eq. 0. Both p c {T) and 
Pab{T) were computed using Montgomery's algorithm. The 
single crystals (A and B) were grown by a flux method based 
on BaO. The ratio a/c for sample B is approximately 8. 



FIG. 2. The resistivity of the crystalline YBa2Cu307 
(Sample C: with volume of unit cell = 390 x 400 x 25 pm) of 
123 phase were obtained by utilizing Eq. [l(j. Both p c (T) and 
Pab{T) were computed using Montgomery's algorithm. The 
ratio a/c for sample C is approximately 16. 
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FIGURE 1 
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Table 1 
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Samples, [Ref.] 


A 


E, 
[J] 


Ei 
[meV] 


[m 2 ] 300 K 
xlO 16 


^Jnp 
[m 2 ] 150 K 
xlO 16 


dp(T) 1 AT 
[l^n-cmK 1 ] 


YBa 2 Cu 3 7 Sample B, [26] 


1.1458 x 10" 8 


89.2£ s 


7.69 


8.022 


5.959 




YBa 2 Cu 3 7 Sample C, [26] 


1.1268 x 10" 8 


94.0^ B 


8.11 


7.895 


5.771 




GaSr 2 (Y . 6 Ca . 4 )Cu 2 O 7 , [27] 


1.3002 x 10" 5 


693k B 


5.98 


8.572 


6.804 


~0 


YSrBaCu 3 7 Sample [O], [13] 


1.1712 x 10 -6 


53.5k B 


AM 


9.036 


7.560 


1.10 


Y . 8 Pr . 2 SrBaCu 3 O 7 Sample [O], [13] 


1.8589 x 10" 6 


76.6k B 


6.61 


8.366 


6.481 


1.60 


YBa 2 Cu 3 O x , [28] 


5.5747 x 10 -6 


95.6k B 


8.25 


7.853 


5.710 


4.25 


YBa 2 Cu 3 7 . x Sample A, [29] 


6.3767 x 10" 6 


92.0k B 


7.94 


7.948 


5.849 


5.00 


Y(Ba . 85 Pb 015 ) 2 Cu 3 O 7 . x Sample B, [29] 


5.4926 x 10 -6 


19M B 


6.88 


8.278 


6.344 


4.90 


Tl(Sr L6 Sm .4)CaCu 2 O 7 . x , [30] 


2.8695 x 10" 5 


59.5k B 


5.13 


8.857 


7.264 


2.60 


Tl(Sr,. 6 Dy .4)CaCu 2 O 7 . x , [30] 


2.2192 x 10" 5 


66.2k B 


5.71 


8.661 


6.946 


2.00 
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